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Abstract 

We have measured the ratio R^ = BR(</> — ► r/7)/BR(0 — ► 777) by looking for 
the radiative decays <f> — > 7/7 and <f> — > 777 into the final states 7r + 7r~77 and 77, 
respectively, in a sample of ~1.4 x 10 9 4> mesons produced at the Frascati <f> factory. 
We obtain = (4.77 ± 0.09 sta t ± 0.19 syst ) x 10" 3 , from which we derive BR(</) -► 
7/7) = (6.20±0.11 s tat ±0.25 sys t) x 10~ 5 . Assuming the rj has zero gluonium content, 
we extract the pseudoscalar mixing angle in the quark-flavor basis, pp = (41.4 ± 
0.3 s tat ± 0.7 sys t ± 0.6th)°- Combining the value of with other constraints, we 
estimate the fractional gluonium content of the rj to be Z 2 = 0.14 ± 0.04 and the 
mixing angle to be tpp = (39.7 ± 0.7)°. 

Key words: e + e~ collisions; radiative (f> decays; pseudoscalar mixing angle 
PACS: 13.65.+i, 14.40. Aq 



1 Introduction 

The value of the 77-77' mixing angle in the pseudoscalar meson nonet, 6p, has 
been discussed extensively over the last 35 years, and is today one of the most 
interesting SU(3)-breaking hadronic parameters to measure [1]. In the context 
of chiral perturbation theory, it has been demonstrated that a description of 
the 77-77' system beyond leading order cannot be achieved in terms of just one 
mixing angle [2]. However, in the flavor basis the two mixing angles are equal, 
apart from terms which violate the OZI-rule; thus working in this basis it is 
still possible to use a single mixing angle, (pp [3]. The ratio of the two 

1 Corresponding author. 
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branching ratios — > 7/7 and — > 777 can be related to the 77-7/' mixing angle 
in the flavor basis [3,4,5,6,7,8] and to the gluonium content of the 77' meson 
[5,9]. 

In this work we present a measurement of obtained using an integrated 
luminosity of 427 pb" 1 collected by KLOE during the years 2001-2002. The 
best measurement available to date was also published by KLOE [10], and 
was obtained from the analysis of the — > 7/7 decay with the 7r + 7r~37 final 
state and the — > 777 decay with the 7r + 7r~7r°7 final state. The earlier KLOE 
measurement was based on an integrated luminosity of about 16 pb _1 , col- 
lected during the year 2000. Previous measurements were performed by the 
SND and CMD-2 collaborations [11,12]. 



2 DA$NE and KLOE 



The Frascati factory DA$NE is an e + e~ collider running at a center-of- 
mass energy y^s = 1.02 GeV, where mesons are produced with the beams 
colliding at a small crossing angle (25 mrad). The KLOE detector consists 
of two main subdetectors: a large cylindrical drift chamber [13] and a sam- 
pling lead/scintillating- fiber electromagnetic calorimeter [14]. A superconduct- 
ing coil surrounding the calorimeter provides a solenoidal field of 0.52 T. 

The drift chamber is 3.3 m in length and has a 2 m radius; it has full stereo 
geometry and operates with a gas mixture of 90% helium/10% isobutane. 
The momentum resolution is cr(p±)/(p±) < 0.4%; the spatial resolution is 
a xy ~ 150 jum and a z ~ 2 mm. Vertices are reconstructed with a spatial 
resolution of ~ 3 mm. 

The calorimeter is divided into a barrel and two endcaps; it surrounds the drift 
chamber and covers 98% of the solid angle. The calorimeter is segmented into 
2440 cells of cross section 4.4 x 4.4 cm 2 in the plane perpendicular to the fibers. 
Each cell is read out at both ends by photomultiplier tubes. Particle arrival 
times and the three-dimensional positions of the energy deposits are obtained 
from the signals collected at the two ends; cells close in time and space are 
grouped into a calorimeter cluster. The cluster energy E is the sum of the cell 
energies, while the cluster time T and position R are energy-weighted averages. 
The probability for a photon to fragment into two or more clusters (splitting) is 
taken into account during event reconstruction by a dedicated procedure. The 
energy resolution is <Je/E = 5.7%/ y^(GeV); the timing resolution is a t = 57 

ps / \jE(GeV) © 100 ps. The KLOE trigger [15] is based on calorimeter and 
chamber information. 
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3 Event selection 



The analysis has been performed using a sample of ~1.4 x 10 9 mesons 
collected in 2001-2002. We search for two categories of events. The first, — > 
7/7 — > 7r + 7r~77, is produced through two different decay chains: 



— > 7/7 with rf — > 7r + 7r 77 and 7/ — > 7r 7r 7r° (1) 

"~ Vl with 77' — > 7T 7r°7/ and 7/ — > 7r + 7r~7r°. (2) 

The second, ^ 777 ^ 77, is produced through 

— > 7/7 with 7/ — ► 7T 7r 7r°. (3) 



The latter channel is used for normalization of the rates and is practically 
background free [16]. Processes with — > KsKl, where the Kl decays near 
the beam interaction point (IP) can mimic the final state with 7r + 7r~77 because 
of the presence of an additional photon due either to machine background or 
split clusters. The most important decay chains of this type are the following: 

-> K S K L with K s -> 7T+7T- and K L -> 7r 7r 7r (4) 
-> i^i^ with X 5 -> 7T°7r and fsT L -> tt+tt^tt (5) 
-> X 5 #l with K s -> 7T+7T-7 and fT L -> 7r 7r 7r (6) 

In particular, the last of these channels has the same final state used to identify 
— > 7/7. Backgrounds from other sources are negligible. 

All events must pass a first-level selection (FLS) [17], consisting of a machine 
background filter and an event selection procedure that assigns events into 
categories. After the FLS, we select 777 (normalization) events by requir- 
ing the presence of seven clusters not associated to tracks. Each cluster must 
have E y > 10 MeV and polar angle 21° < # 7 < 159° with respect to the beam 
direction, in order to exclude the focusing quadrupoles. Each cluster time must 
be compatible with the hypothesis of a photon coming from the IP, so that 
the requirement to be satisfied is |T 7 — Ry/c\ < 5<7 t , where T 7 and i? 7 are 
the cluster time and position and a t is the time resolution. We further require 
that no tracks come from the interaction region. To select — > 7/7 events, we 
similarly require seven prompt photons; in addition, we require the presence 
of a vertex formed by two tracks of opposite charge within a cylindrical region 
\J x ltx + 2/vtx < 4 cm and \z vtx \ < 8 cm around the IP. After applying these 
selections, we perform a kinematic fit requiring energy-momentum conserva- 
tion, and times and path lengths to be consistent with the speed of light for 
photon candidates. We perform an inclusive measurement of the two different 
decay chains contributing to the — > 7/7 process. 

The background is significantly reduced by the event-classification cuts in the 
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FLS, which efficiently identify K$ decays. By selecting events classified only 
as pure radiative decays, the K S K L backgrounds are largely reduced. After 
this cut, background from processes (4) and (6) above is practically negligible, 
while to further reduce background from process (5), we require E 1 > 20 MeV 
for each of the seven clusters, where the improved estimates of the cluster 
energies from the kinematic fit are used to evaluate this cut. 

To estimate the absolute number of residual background events, the number of 
<f> mesons in the sample has been evaluated. We obtain Nj, = (1.37±0.02) x 10 9 , 
where the uncertainty comes from the estimate of the integrated luminosity 
and the experimentally-determined value of the leptonic decay width of the 
(f). The luminosity is known with a total error of 0.6% [18], while the leptonic 
width has been measured by KLOE with a 1.7% uncertainty [19]. 

For each background channel, Table 1 lists the branching ratio multiplied by 
the overall selection efficiency as determined by Monte Carlo (MC) simula- 
tion. The branching ratio for process (6) has been evaluated using BR,(Ks — > 
7T+7T-7) = (4.87 ± 0.11) x 10~ 3 , for £ 7 > 20 MeV/c [20], the same cut used 
for our analysis. 

Residual background due to multiphoton final states with a photon conversion 
and production of a charged-particle vertex has been reduced to a negligible 
level by a cut on the opening angle of the two tracks. Requiring cos 6 l vr + vr - < 
0.84 [16] leads to a small effect on the signal efficiency and gives rise to a 
negligible contribution to the systematic error. 



Table 1 

Branching ratios times selection efficiencies BR(^ — > KsKl) ■ BR(/fs — ► 7T7t(7)) • 
BR(A'i — > mm) ■ e(MC) and event yield for <f> —* K$Kl backgrounds 



K S /K L - 


^+^-/vr 7r 7r° 




7T + 7T 7/7T 7r 7T 


BRxe(MC) 


(44 ± 3) • 10~ 9 


(112 ±4) • 10~ 9 


(95 ± 3) • 10" 9 


Background events 


60 ±4 


153 ±6 


130 ±5 



After all selection cuts, we find 3750 events with an expected background 
of ^background = 343 ± 43. The background uncertainty comes from the un- 
certainties on the number of mesons, the selection efficiency (±15), and 
the evaluation of the machine background level (±28) [16]. We treat the back- 
ground uncertainty as a systematic contribution to the error on N v >^. The final 
number of events from processes (1) and (2), after background subtraction, is 
Nrfj = 3407 ± 61 sta t ± 43 syst . The event- select ion efficiencjd] after FLS cuts is 

2 We perform an inclusive measurement of two different decay chains, (1) and (2), 
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evaluated by MC simulation to be e v > = (23.45 ± 0.16)%. The sample of MC 
signal events is five times larger than the data sample and has been produced 
by carefully simulating the run conditions during 2001-2002 data taking. Cor- 
rections for small data-MC discrepancies have been estimated directly from 
control data samples and applied. The FLS efficiency has been measured from 
data using a small subsample of minimum-bias events reconstructed without 
the FLS filters; the result is e^ LS = (97 ± 1)%. 

Figure 1, left, shows the data-MC comparison of the photon-energy spectrum 
after all cuts. For 7/7 events, the agreement is good. However, it is not possible 
to identify the recoil photon since its spectrum overlaps with the spectrum for 
photons from 7/ decay (see Fig. 1, right). Therefore, when building the 7T + 7T~67 
invariant mass, we have a combinatorial background due to incorrect photon 
association. If we subtract the MC estimate of the physical and combinatorial 
background from data, we obtain a peak at the expected value of the 7/ mass, 
as seen in Fig. 2. 
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Fig. 1. <p — > 7/7. Energy distribution for all seven prompt neutral clusters. Left: 
Data-MC comparison after all cuts. Dots show data; histograms show MC distribu- 
tions for each of the three background processes, (6), (5), and (4), and for the sum 
of signal and background contributions. Right: MC distributions for the recoil 7 
emitted in the 4> decay, the six 7's from the 7/ decay, and the overall distribution. 

In contrast to the case of — > 7/7 decay, in — > 777 decay the recoil photon 
is the most energetic photon in the event. As a result, the recoil photons give 
rise to a distinct peak in the photon-energy spectrum. To select process (3), 
we additionally require the identification of the recoil photon with 320 MeV < 
Ey < 400 MeV. This cut is fully efficient for signal events and allows us to build 

that contribute to the <p — > 7/7 process. The event-selection efficiency is obtained 
from the branching-ratio weighted average: e v ' = (e^BR^ + £( 2 )BR( 2 ))/(BR( 1 ) + 

BR (2)) 
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Fig. 2. <j> — > 7/7. Left: Distribution of M{k + it~Q~{) for all possible combinations. 
Dots are data; the histogram represents the MC sum of the three background con- 
tributions, (6), (5) and (4), the combinatorial background, and the contribution 
from correctly paired signal events. Right: Data-MC comparison of M(ir + n~Q~{) 
distributions in the area about the peak, after bin-by-bin subtraction. 

the 67 invariant mass, as seen from Fig. 3. In all, we select N V1 = 1 665 000 
events. The event selection efficiency after FLS is evaluated from MC; its 
value is e v = (33.66 ± 0.01)%. The FLS efficiency measured using data is 
e? LS = (97.88 ±0.10)%. 
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Fig. 3. <p — ► r/7. Left: Data-MC comparison of distributions. Right: Data-MC 
comparison of M(6j) distributions. In both panels, histogram shows MC; dots are 
data. 
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4 Results 



We evaluate the ratio of the two branching fractions R<j, = BR(0 — > 7/7) /BR(0 - 
777) using the relation 

p Nyh £ v r< BR^ is / 7 \ 

where BR^ is the branching fraction BR(7/ — > 7r 7r 7r ) and BR^/ 7 is BR(t/ — > 
tt+tt-^-BR^ -> ttW) + BR(t/ -> ttVV) • BR(?7 -> tt+tt-tt ); e v and ^ are 
the overall selection efficiencies respectively for — > 7/7 and for — > 777; and 
Cfls — e ^ LS / e % LS is the ratio of the FLS efficiencies. The factor K p accounts 
for interference between the amplitudes A(<f> — > r)(r)')~f) and A(p — > ^(rj')^), 
and has been evaluated in the Born approximation in the manner of Ref. [10]. 
This correction depends on y/s, for which a mean value over run conditions 
has been evaluated; the standard deviation is taken as a contribution to the 
systematic uncertainty. The final value K p = 0.95 ± 0.01 takes into account 
also radiative corrections to the cross section [19]. 

The systematic uncertainty on is dominated by the uncertainty on the in- 
termediate branching ratios of the 7/. The evaluation of the systematic uncer- 
tainty on the ratio of analysis efficiencies is obtained by studying the samples 
and the stability of the results with respect to (i) the first level selection, (ii) 
the tracking and vertex efficiencies, and (iii) the cuts applied to select the 
final sample. The contribution to the uncertainty due to FLS comes from the 
uncertainty on the efficiency measured from data using a small subsample of 
minimum-bias events reconstructed from raw data without the FLS filters. 
The error due to the tracking and vertex efficiencies is estimated from the 
data-MC discrepancy observed for a control sample of — > it + it~-k events 
with charged-pion momenta in the range covered by our analysis. A sizeable 
contribution to the uncertainty due to the final selection cuts comes from the 
inclusiveness of the cut on the \ 2 probability from the kinematic fit, since the 
xV^Vdof distribution has a pronounced tail. We estimate this effect by moving 
the cut and evaluating the maximum variation of the ratio A(N/e)/(N/e) and 
its effect on R^. The systematic error on N^^/N^ comes from the background 
subtraction. All the contributions are summarized in Table 2. 

The final result is: 

= (4.77 ± 0.09 stat ± 0. 19 syst ) x 10" 3 (8) 

Using the value of BR(0 -> 777) = (1.301 ± 0.024)% from the PDG [21] we 
find: 

BR(0 -> 7/7) = (6.20 ± 0.1W ± 0.25 syst ) x lO" 5 (9) 
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Table 2 



Contributions to the systematic error on 



Quantity with systematic source 


Systematic error 


Source 


N rh /N vl 


1.3% 


Background 


Cfls 


1.0% 


Preselection 




1.0% 


Tracking efficiency 


Sjj/ Erf 


1.0% 


Vertex efficiency 




1.5% 


Prob(x 2 ) 


BR^/BR^/^ 


2.8% 






1.0% 




Total 


4% 





The value of can be related to the pseudoscalar mixing angle. In the quark- 
flavor basis, the 77-77' system can be parameterized in terms of just one angle: 

\rf) = coscpp \qq) + siny?p \ss) , 
= — sinipp \qq) + cos ip P \ss) . 

where \qq) = (1/^/2) uu + dctj. Using the approach of Ref. [7,8], where SU(3)- 
breaking is taken into account via the constituent-quark mass ratio m s /rh, R<f, 
is given by the following expression: 

BR(0^r/7) 2 / m s C NS tan ip v \ 2 (p v ,\ 3 

R <P = "rJrJT^ 7 = cot V?p 1 - —r, — — (10) 

BR(0 -> 777) \ m C s sin 2tp P J \p v J 

where tpy = 3.4° is the mixing angle for vector mesons, Pm) is the 77(77') 
momentum in the <fi center-of-mass, and the two parameters C^s an d C$ (see 
Table 3) represent the effect of the OZI-rule, which reduces the vector and 
pseudoscalar wave-function overlap [8]. From Eq. (10) we obtain the following 
result: 

V? P = (41.4 ± 0.3 stat ± 0.7 syst ± 0.6th)°. (H) 
The theoretical uncertainty on the mixing angle has been evaluated from the 
maximum variation induced by the spread of the values for m s /m, Cats, and 
Cs- In the traditional approach, 77-77' mixing is parameterized in the octet- 
singlet basis; in this basis the value of the mixing angle becomes: 9 P = ip P — 
arctan ^2 = (-13.3 ± 0.3 stat ± 0.7 syst ± 0.6 th )°. 

In QCD, gluons may form a bound state, called gluonium, that can mix with 
neutral mesons. While the 77 is well understood as an SU(3)-flavor octet meson 
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with a small singlet admixture, the rj' is a good candidate to have a sizeable 
gluonium content. If we allow for non-zero rj' gluonium content, we have the 
following parameterization: 

\rj') = X v > \qq) + Y n > \ss) + |gluons) 

where parameterizes mixing with gluonium. Normalization implies X 2 , + 
Y 2 + Z 2 =1 with 

X v > = cos0Gsiny9p 

Y r f = cos (f) G cos (f P (12) 
Z v > = sin^G 

where (fie is the mixing angle for the gluonium contribution. Non-zero rj' glu- 
onium content would imply 

X 2 V , + Y 2 <1 (13) 
and Eq. (10) would have to be rewritten 

d ^ 2, (-, m s C NS t&nip v \ 2 (p v ,\ 3 

i? = cot y? P cos G 1 - r-^ — — (14) 

\ m C s sm2(f P J \p v J 

We may use SU(3) relations as proposed in Refs. [5,8,9] to further constrain 
X v > and Y v >. With the parameterization of Eq. (12), these relations may be 
written as: 

F(rj' — > 77)/r(7r° — > 77) = - (5 cos 4>g sin (f P + V2^- cos 4>g cos (fp) 2 

9 \m n J ' ts 



C N c „ (ml,- m 2 m u 



W ' \ /TV \ ^NS o / "V "V " l u \ 2 i • 2 

r(?7 -> fn)/T(u -> tt u 7) = 3 £ cos G sm ip P 

cos 991/ \ m~ - m% m v > J 

(2 2 \ ^ 

m v ,-m rn^A ^ ^ 

- Try / 

777- 

+ • tan^y • cos0g cosifp} 2 

m 

Figure 4 shows the above constraints graphically in the (X v r,Y v r) plane. The 
constraint from the KLOE measurement of is shown in the hypothesis 
of no rj' gluonium content. To allow for gluonium, we have minimized a \ 2 
function with cos 2 0g and cos 2 <^p as free parameters by imposing the above 
constraints and including in the error matrix the uncertainties on the other 
parameters as shown in Table 3. The values for f q and f s are from Ref. [3]; all 
other values are from [8] . We assume no correlations between the parameters. 

The solution in the hypothesis of no gluonium content, i.e. cos 2 G = 1, yields 
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Fig. 4. Bounds on Xj, Y v >. (1): T(rj' -» pj)/T(oj -> vr°7); (2): 

r(V - 77)/r(^° -> 77); kloe ^ r(0 -> r/ 7 )/r(0 - r? 7 ); (3): 
r(?/ -> w 7 )/r(o; -► 7r°7). 



Table 3 



Fixed parameters used 


in the fit for 


cos 2 99(j and 


COS 2 ifp 




Parameter / 9 


/, 


Cns 


Cs 


rn s 
rh 


Value 1 ± 0.01 


1.4 ±0.014 


0.91 ±0.05 


0.89 ± 0.07 


1.24 ±0.07 



ifp = (41.5*%%)°. The value of xV-^dof is 11-34/3, corresponding to a proba- 
bility of 0.01, and suggests a possible non-zero value for Z 2 ,. 

The solution allowing for gluonium is cos 2 <Pg — 0.86 ± 0.04 and cos 2 <y9p = 
0.592 ± 0.012, from which <p P = (39.7 ± 0.7)° and Z 2 ql = 0.14 ± 0.04, which 
means \ipa\ = (22 ± 3)°. In this case, x 2 /N Aoi = 1.42/2, for a probability of 
0.49. When the constraint from rj' — > wy decay is relaxed, similar results are 
obtained: <p P = (39.8 ± 0.8)° and Z 2 = 0.13 ± 0.04, while the \ 2 probability 
improves to 0.97. Figure 5 shows the four constraints in the (Z 2 ,,{p P ) plane, 
together with the allowed region for the solution with gluonium. 

In the octet-singlet basis, the value of the rj-rj' mixing angle fp = (39.7 ± 
0.7)° becomes dp = (—15.0 ± 0.7)°, which is in agreement with the value 
9p = (—15.9 ± 1.2)° measured by BES-II [22] using the branching fractions 
BR(J/ijj — > 7/7) and BR(J/ip — > 777). These results for ip P are consistent 
with those coming from other processes as reviewed in Ref. [3], which quotes 
an average value <fp = (39.3 ± 1.0)°. 
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Fig. 5. The four constraints in the (Z 2 ,(fp) plane: y\ = T(r]' — ► 77)/r(7r° — > 77), 
V2 = r(V p7)/r(w 7T°7), y 3 = i^, and y 4 = r(r/ w 7 )/T(u; 7T°7). The 
allowed region for the solution with gluonium is shown. 

5 Conclusion 

Analyzing about 1.4x 10 9 mesons collected by KLOE at the DA$NE collider, 
we have obtained a precise new measurement of the ratio = BR(0 — > 
r/7)/BR(0 -> 777) = (4.77±0.09 st at±0.19 syst ) x 10~ 3 . From this measurement, 
we obtain the most precise value of the pseudoscalar mixing angle in the flavor 
basis: ipp = (41.4 ± 0.3 sta t ± 0.7 syst ± 0.6 t h)°- Allowing for gluonium in the 77', 
a fit to our result together with other measurements yields ipp = (39.7 ±0.7)° 
and a 7/ gluonium content of Z 2 , = 0.14 ± 0.04. 
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